Coulomb dissociation, a tool for nuclear astrophysics 
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O . 1. Introduction 
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^ I Cross sections for nuclear reactions of astrophysical interest can be very small, too small 
O ' to be directly accessible experimentally. Sometimes indirect methods can be useful to 
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determine astrophysical S-factors. The ANC (Asymptotic Normalization Constant) 
method [1], the Trojan horse method [2], [3l HJ [5] and the Coulomb dissociation method 
[6] were discussed at the present conference. A more recent minireview is given in 



^ ' [7]. In all these methods, some theoretical interpretation of the experimental results is 

necessary in order to arrive at an astrophysical S-f actor. 

In the Coulomb dissociation method a fast projectile passes through the Coulomb 
field of a heavy nucleus. The time-varying electromagnetic field experienced by the 
projectile acts like a spectrum of equivalent photons which can dissociate the projectile. 
The equivalent photon spectrum can be computed from the kinematics of the process. 
In this way one can determine photodissociation cross-sections, essentially based on 
quantum electrodynamics. Radiative capture cross sections leading to the ground state 
of the final nucleus (the projectile) can be deduced by detailed balance. The projectile 
can be an exotic nucleus. Thus one has the unique possibility to study the interaction of 
unstable nuclei with photons [HI [9]. In particle physics the Primakoff effect [10] has been 
used for some time to study the interaction of photons, pions, A's, etc. with photons. 
Another method to study the interaction of photons with exotic nuclei are e-A colliders 
Coulomb dissociation is the only method for many years to come. 
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The theoretical description of intermediate energy Coulomb excitation and 
dissociation has been developed over the past decades, it is reviewed in ^2]. The main 
characteristics are given in Sect. 2. We propose as a 'homework problem' a quasirealistic 
simple model, which can be solved analytically. It could also serve as a benchmark for 
tests of more sophisticated or involved models like the Continuum Discretized Coupled 
Channels model (CDCC). A major result of radioactive beam physics is the discovery of 
low-lying electric dipole strength in neutron-rich nuclei, see, e.g., [13]. This discovery was 
made possible by the 'working horse' Coulomb dissociation. This is discussed in Sect. 3. 
Low lying dipole strength is directly related to the halo structure of these nuclei. The 
effective-range approach to low lying El strength and simple scaling laws are discussed. 
In some astrophysical scenarios for the r-process it is vital to know the low lying electric 
dipole strength. With the future radioactive beam facilities it will be possible to access 
these questions. 



2. Theory of intermediate energy Coulomb excitation and dissociation 

One of the basic parameters in Coulomb excitation is the ratio of the collision time to 
the nuclear excitation time, the so-called adiabaticity parameter 

( = - (1) 

where ut is the excitation energy, b is the impact parameter and v the beam velocity. 
The corresponding Lorentz factor 7 is typically not much larger than one. For ^ <^ 1 the 
process is sudden, and excitation is possible; for ^ 1 the system follows adiabatically 
the time varying field and the excitation probability tends to zero. The strength of the 
excitation is measured by the strength parameter 
(,) ^ ZreiflMjEm 

^ hvb^ ■ 

The projectile, characterized by the electromagnetic matrixelements {f\M{EX)\i) is 
excited on its passage through the Coulomb field of the (heavy) target nucleus with 
charge Zt- The parameter x^'^'' can be regarded as the number of exchanged photons. 
In electromagnetic excitation it is a good approximation that the nuclei do not penetrate 
each other. In this case, the interaction is encoded in the electromagnetic matrixelement 
{f\M{EX)\i) between the relevant states. The Coulomb parameter rj is the monopole 
strength parameter, i.e. A = in ([2]), and the multipole matrixelement is replaced by 
the charge Ze of the nucleus. For higher beam velocities higher order effects tend to 
be small. For not too light nuclei 1] is still ^ 1, and the semiclassical description is 
appropriate. 

2.1. A quasirealistic and analytically solvable model of Coulomb excitation of neutron 
halo nuclei 

An archetype of a halo nucleus is ^^Be with a ^°Be core and a single halo neutron 
in the 2si/2-state. There is a strong El transition to the l/2~ bound state, the only 
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bound excited state of the system. This dipole transition was studied by Coulomb 
excitation at GANIL, RIKEN and MSU [H [IS [ISl [H] . In Ref. higher order effects 
are studied in the 'XCDCC approach. We take this as an opportunity to revisit the 
theoretical study of higher order effects in intermediate energy Coulomb excitation. In 
[18j electromagnetic excitation of ^^Be is studied in the sudden limit of the semiclassical 
method. Higher order effects are treated to all orders. In [19] an analytically solvable 
model for higher order electromagnetic excitation effects of neutron halo nuclei was 
presented. In that work, there was only the s-wave bound state. Now we consider the 
case where there is, in addition, a p-wave bound state, as it is the case in ^^Be. In the 
sudden approximation the dipole excitation amplitude is given by 

flsuddcn = (/I exp (-igcoul ■ ^) (3) 

where gcoui = ^^^^^'^x- The impact parameter b is chosen to be in the x-direction. 
(For the actual calculation in polar coordinates it is convenient to change this to the 
z-direction.) The dipole approximation is quite well fulfilled, the dipole effective charge 
■^cff — m^+m much larger than the corresponding quadrupole charge. The neutron 
and core mass are denoted by m„ and iric respectively, the charge of the core is given 
by Zc- The sudden approximation is applicable for ^ <^ 1. Even for the comparatively 
low GANIL energies of about 40 MeV/nucleon this is reasonably well fulfilled. The 
most important intermediate states are expected to be in the low energy continuum, 
where the dipole strength has a peak, at around 1 MeV excitation energy. The sudden 
approximation has the advantage that intermediate states are treated by closure, thus 
one only needs a model for the initial and final states, and not for all the intermediate 
states. In lowest order in gcoui the first order dipole approximation is obtained. It is 
shown in [18] that third order El excitation is more important than second order E1-E2 
excitation. The matrixelements are dominated by the exterior contributions. In a pure 
single particle model the radial wave functions of the 2s 1/2 and lpi/2 states are given by 

/o(r) = CQq^rhaiiq^r) (4) 

and 

/i(r) = Ciqirhiiiqir) (5) 

Both states are halo states and the normalization constants are given by Cq = -^/Sgo 

and Ci = ^2qiR/3 in the halo hmit respectively [2Cy. The bound state parameters qi 

^22 

{i = 0, 1) are related to the binding energies by Ei = where is the reduced mass of 
the core-neutron system. We have Eq = 504 keV and Ei = 184 keV. The radius of the 
core is denoted by R. With these model asumptions we can calculate the B(E1) value 
for the 1/2+ 1/2 ^-transition as well as the higher order effects in electromagnetic 
excitation. Whereas in [T7] quite sophisticated models are used, our approach is simple, 
transparent and at the same time close to reality. We propose this to be a model study 
and leave the spectroscopic factors equal to one. (They could be adjusted, which would 
result in a quasi- reahstic description of the ^^Be system for our purpose.) The XCDCC 
calculations are quite involved, with many parameters. It would be very useful to check 
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the method by comparing to a simple case, such as this one, where analytical results are 
possible. In order to avoid nuclear effects a sharp cutoff at a minimum impact parameter 
femin can be introduced. 



The B{E1) value is given by B{E1) = (z^\\eY |i?oi| V(47r), where Rqi is the radial 



, 2 

dipole integral. The integrals are elementary. Using /q°° r"e~"^dr = n\/a"'~^^ we find 
270^/^(70 + 271) 

where 70 = qoR and 71 = qiR. We note that we extended the radial integral over 
the exterior wave function from R to zero. For i? — the radial dipole integral goes 
to zero because the normalization of the p-wave function tends to zero in this limit. 
Thus R must be kept finite, say R = 2.78 fm [21], this value determines the asymptotic 
normalization of the p-wave bound state. We find B{E1) = 0.193 e^fm^, to be compared 
to the value of B{E1) = 0.105(12) e^fm^ obtained from an analysis of the GANIL data, 
see [T7j, consistent with other Coulomb dissociation experiments at RIKEN and MSU 
and the value obtained by the Doppler shift attenuation method [22] . 

We expand the excitation amplitude ([3]) in terms of the dimensionless strength 
parameter y = qcon\/{.% + Q'l)- (Cf. ([2]), we take e/(go + ?i) as a convenient measure 
for the order of magnitude of the dipole matrix-element.) The excitation probability 
is given by P{h) = |asuddenP- The lowest order term is proportional to y"^. The most 
important higher order contribution comes from the third order in gcoui- It can be 
calculated analytically. Its interference with the lowest order term leads to the next 
term in the expansion in y, of the order of y'^. We have P{b) = Plq + Pnlo + •••• The 
lowest order term is given by 

_ , 470(70 + 271)' _ ^ /,2 (J^ 

27(7o + 7i) 
The next term is found to be 

p ^ 4 870^^(70 + 27i)(7o + 471) _ ^ /,4 (c^^ 
P^^o = -y 45(^^ + ^^)2 = -^4/6 (8) 

Total cross sections are obtained by integration over the impact parameter, starting 
from a minimum impact parameter 6min- The sudden approximation fails for large 
impact parameters, and an adiabatic cut-off femax = 'jv/u has to be introduced for the 
lowest order result. We put lu = 320 keV, the energy of the 1/2" state in ^^Be. (For the 
higher order terms this is not necessary, the convergence in b is fast enough.) We get 

aLO = 2nC2 In ^ (9) 



and 



7rC4 

(^NLO = -To • (^0) 



,2 
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We note that the strength parameter y is proportional to 1/v, i.e. the leading order 
term decreases like l/E, the next-to-leading order term like l/E"^, where E is the beam 
energy. We think that this analytical model could serve as a benchmark for tests of more 
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involved reaction models. We hope to publish a more detailed account of the present 
approach in the future. 



3. Low-lying electric dipole strength in neutron rich nuclei 

An effective-range approach to low lying El-strength for one-neutron halo nuclei was 
developed in [201 121] . There is a small parameter 



7^qR=^- (11) 

In lowest order, the dipole strength is independent of 7. The S(i?l)-strength function 
is proportional to the shape function S^^ and scales with the parameter = E/Eund = 
IQqi where E is the cm. energy in the continuum. For s-p transitions it is given by 

m 

= {ilx^Y " ""'^'^^ + ^"^'^^^ + •••] • ^^^^ 

This remarkably simple result can be apphed to deuteron photodisintegration, and si/2- 
neutron halo nuclei like ^^Be, ^^C, ^^C,... The interaction of the final state p-wave 
neutron with the core can usually be neglected. Thus low lying strength due to transition 
to a structureless continuum is found. It may look like a resonance, but it has nothing 
to do with a resonance. This was recognized long time ago [221 121] • We quote from a 
recent review of low lying dipole strength [26]: '...the onset of dipole strength in the 
low-energy region is caused by nonresonant independent single-particle excitations of 
the last bound neutrons'. In general there are characteristic effects of the core-neutron 
interaction in the continuum state, usually more pronounced for states with If = li — 1. 
For the s-p transitions this term is proportional to 7^, which is quite small for a halo 
nucleus. For the low energies relevant here this interaction can be parametrized in 
terms of the scattering length. An interesting effect of this type was found by analysing 
the high precision data of ^^Be Coulomb dissociation [27]. A large scattering length 
O'iZi^'^ = 456 fm'^ was found [21]. It is due to the pi/2-subthreshold state. A treatment 
of two-neutron halo nuclei in the effective range method for low lying strength of halo 
nuclei is given in [25]. 



4. Conclusion and Outlook 



Electromagnetic excitation is a powerful tool to investigate the interaction of (quasireal) 
photons with unstable nuclei. It will continue to play a prominent role at the future 
radioactive beam facilities. A good theoretical understanding of the process and its 
interplay with nuclear excitation is mandatory, see, e.g., |6|. In the future rp-process 
nuclei will come into focus. The possibility of 2p-capture is also discussed. It will never 
be possible to study this process in the laboratory. However, the time-reversed process 
of Coulomb dissociation with two protons in the final state is well within reach. An 
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example is Coulomb dissociation of ^^Ne where the soft dipole mode in this proton- 
rich nucleus is discussed. By a suitable Coulomb dissociation experiment valuable 
information on the 2p-capture cross section on ^^O at astrophysical conditions approriate 
for explosive burning in novae and X-ray bursts may be obtained. 

At future radioactive beam facilities r-process nuclei will become available. In 
certain scenarios (see, e.g., it will be important to know the low lying El strength, 
which will decisively influence the r-process abundances. 
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